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FENG DAI AND YUAN XU 



Abstract. This is Chapter 1 of the book Approximation Theory and Har- 
monic Analysis on Spheres and Balls by the authors. It provides a self- 
contained introduction to spherical harmonics. The book will be published 
as a title in Springer Monographs in Mathematics by Springer in 2013. The 
table of contents of the book is attached at the end of this file. 



In this chapter we introduce spherical harmonics and study their properties. 
Most of the material of this chapter, except the last section, is classical. We strive 
for a succinct account of the theory of spherical harmonics. After a standard treat- 
ment of the space of spherical harmonics and orthogonal bases in the first section, 
the orthogonal projection operator and reproducing kernels, also known as zonal 
harmonics, are developed in greater details in the second section, because of their 
central role in harmonic analysis and approximation theory. As an application of the 
addition formula, it is shown in the third section that there exist bases of spherical 
harmonics consisting of entirely zonal harmonics. The Laplace-Beltrami operator 
is discussed in the fourth section, where an elementary and self-contained approach 
is adopted. Spherical coordinates and an explicit orthonormal basis of spherical 
harmonics in these coordinates are presented the fifth section. These formulas in 
two and three variables are collected in the sixth section for easy reference, since 
they are most often used in applications. The connection to group representation 
is treated briefly in the seventh section. The last section deals with derivatives and 
integrals on the sphere. With the introduction of angular derivatives that are first- 
order differential operators acting on the large circles of intersections of the sphere 
and the coordinate planes, it is shown that the Laplace-Beltrami operator can be 
decomposed into second-order angular derivatives. These derivative operators will 
play an important role in approximation theory on the sphere. They are used to 
derive several integral formulas on the sphere. 

1. Space of spherical harmonics and orthogonal bases 

We begin by introducing some notation that will be used throughout this book. 
For X G M.'^, we write x — {xi, . . . , Xd)- The inner product oi x,y G K'' is denoted 
by {x,y) := X]i=i ^i?/* ^^^ the norm of x is denoted by ||a::|| := •^/(a;, x). Let Np 
denote the set of nonnegative integers. For a = (ai, . . . , ad) & Nq, a monomial x" 
is a product x" = x"^ . . . x'^'' , which has degree \a\ = ai + . . . + ad- 
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A homogeneous polynomial P of degree n is a linear combination of monomials of 
degree n, that is, P{x) — 'J2\a\=n Cqs;", where Cq are either real or complex numbers. 
A polynomial of (total) degree at most n is of the form P{x) — J2\a\<n Cax""- Let 
Vn denote the space of real homogeneous polynomials of degree n and let IIJ^ denote 
the space of real polynomials of degree at most n. Counting the cardinalities of 
{a e Ng : \a\ == n} and {a G Ng : \a\ < n} shows that 

dim-p^= and dimOf = 



n J \ n 

Let di denote the partial derivative in the i-th variable and A the Laplacian 
operator 



A:=9? + ... + 9^ 



Definition 1.1. For n = 0,1,2,... let Tif^ be the linear space of real harmonic 
polynomials, homogeneous of degree n, on M."^ , that is, 

nt:={Per^:AP = 0}. 

Spherical harmonics are the restrictions of elements in H^ to the unit sphere. If 
Y e H^, theny(a;) = ||a;||"y(a:;') where a; = ||2;||a;' and x' G S'^^^. Strictly speaking, 
one should make a distinction between "H^ and its restriction to the sphere. We 
will, however, also call "HJ^ the space of spherical harmonics. When it is necessary 
to emphasize the restriction to the sphere, we shall use the notation Hn\sd-i. In 
the same vein, we shall define P„(S''"^) := P^|s<i-i and n„(§''-i) := n^|sd-i. 

Spherical harmonics of different degrees are orthogonal with respect to 

(1.1) {f^g)^,_,:=±f f{x)g{x)d<j{x) 

^d Jgd-l 

where dcr is the surface area measure and ujd denotes the surface area of S'^^^, 



(1.2) UJd:= / da = 



d-1 



T{d/2y 

Theorem 1.2. IfYn G 71^, Ym G 'H'^i, and n^ m, then (y„,y„j)gd_i — 0. 

Proof. Let ^ denote the normal derivative. Since Yn is homogeneous, Yn{x) = 
r"r„(x'), where x = rx' and x' e S'^-i, so that ^(x') = nY,,{x') for x' G '&^-^ 
and n > 0. By Green's identity, 

(n-m) \ YnYmda ^ / [Ym^^ - Y^^^jda 
Jsd~i 7sd-iV or or J 

V„,Ar„-y„Ay„)dx = o, 

since Ar„ = and Ay„i = 0. D 

Theorem 1.3. For n = 0, 1, 2, . . ., there is a decomposition ofV^, 

(1-3) v'^= WxfmU^. 

0<j<n/2 

In other words, for each P G Vf^, there is a unique decomposition 

(1.4) P(.t)- ^ \\x\\^^Pn-2j{x) with Pn-2j(^K-2j- 

0<j<n/2 



Proof. The proof uses induction. Evidently Vq = Hq and Vf = Tif. Since 



AT'^ C Vn-2J diuiHfi > dimT'^ — dim7'^_2- Suppose the statement holds for 



771 = 0, 1, ... , n— 1. Then ||x|p7^^„2 is ^ subspace of P^ and it is isomorphic to 'P^_2- 
By the induction hypothesis, \\x\\^V^_2 = 0o<j<„/2-i bll^^^^'^n-2-2r Hence, by 
the previous theorem, 'H^ is orthogonal to ||a:;|P'P^_2j so that dimH^ + dimP^.g < 
dim-p^. Consequently, V^ = Hi S) \\x\\^Vn-2- ^ 

Corollary 1.4. For n ^ 0,1,2,. . ., 

(1.5) dim-H^ = dim-p,t - dim-p^_2 = 

where it is agreed that dim7'„_2 ^ f^^ n = 0,1. 

Corollary 1.5. ForneN, n„(§''-i) = 'P„(§''"^) ffi Vn-iiS'^^^) and 

(1.6) dimn„(§'^-i) =dim7';J + dimnti = (''^^"^) + C't-i^ 

Proof. By Theorem 1.3, n„(§''~^) can be written as a direct sum of Tif for < 
k < n, which gives the stated decomposition by (1.3). Moreover, 



ri ,. „^ ,. .„w iii + d—l\ in + d — 3 

^ n 



dim n„ (§^- 1 ) = ^ dim -H^ = ^ (dim V^ - dim V^ 



-'k-2) 
fe=0 *:=0 



by (1.5), which simplifies to (1.6). D 

The orthogonality and homogeneity define spherical harmonics. 

Proposition 1.6. If P is a homogeneous polynomial of degree n and P is orthog- 
onal to all polynomials of degree less than n with respect to (•, •)gd-i, then P G H^. 

Proof. Since P € V^, P can be expressed as in (1.4). The orthogonality then shows 
that P = P„e-Ht □ 

Let 0{d) denote the orthogonal group, the group of d x d orthogonal matrices, 
and let SO{d) = {g G 0{d) : det g — l}he the special orthogonal group. A rotation 
in W^ is determined by an element in SO{d). 

Theorem 1.7. The space H^ is invariant under the action f{x) i— >■ f{Qx), Q E 
0{d). Moreover, if {Ya} is an orthonormal basis ofT-V^, then so is {i^Q:(Q{'})}- 

Proof. Since A is invariant under the rotation group 0{d) (writing A = V • V and 
changing variables), if F G Ut and Q e 0{d) then Y{Qx) e Ui. That {Ya{Qx)} 
is an orthonormal basis of T-V^ whenever {Ya{x)} is follows from 

— / Yo,{Qx)Yf,{Qx)da{x) - — / Y^{x)Yfi{x)da{x) = S^.p, 

which holds under a change of variables since da is invariant under 0{d). D 

Besides {f,g)^d-i, another useful inner product can be defined on Vjl through 
the action of differentiation. For a € NJ^, let d" := d"^ ■■■d^''- Let (a)„ :== 
a{a + 1) • • • (a + n — 1) be the Pochhammer symbol. 

Theorem 1.8. For p,q e "P^, define a bilinear form 

(1.7) {p,q)^:^p{d)q, 

where p{d) is the differential operator defined by replacing x" in p{x) by d" . Then 



4 FENG DAI AND YUAN XU 

1. (p,q)g is an inner product on V^; 

2. the reproducing kernel of this inner product is kn{x,y) :~ {x,y) /n!; that is, 

{kn{x,-),p)g^p{x), VpeVi; 

3. forpeV^ andqeni, 

Proof. Let p,q E V^ be given by p{x) = E|a|=« a^^;" and q(x) = E|a|=«^"2;", 
where aa,ba G M. Then, 

(1.8) {p,q)d^ X! "^"^^ X! ^'9^'^= X! Q^'""^", 

|a|— n |/3|— n |a|— n 

which imphes, in particular, that {p,p)g > for p ^ 0. It fohows then that (•, ■)g 
is an inner product on V^. By the multinomial formula, for qa{x) = x", \a\ — n, 

{k^{x,-),q.,)g = - J2 Q^^^^^^^^ala;), 

which shows that kn{x,y) is the reproducing kernel with respect to {•,-)g- 
We now prove item 3. Integrating by parts shows that 



a,/(x)5(x)e-ll^ll /'dx = - / fixMgix) - x,g(x))e-ll-ll /'dx. 
Since p{d)q is a constant, using this integration by parts repeatedly shows that 



(27r)'i/2 
1 



Kx)(p(a;)+s(x))e-ll"ll'/2d^ 



(27r)''/2 

where s E Il'^_^. Since q G "H^ and p G T'^'J, switching to polar integral and using 
the orthogonality of H"^, we obtain 

-I /'OO /• 

(P' '?)a - 7:71^ / r^^+^-ie-V^dr / q{x')p{x')d<j{x'). 
Evaluating the integral in r and simplifying by (1.2) concludes the proof. D 

A large number of spherical harmonic polynomials can be defined explicitly 
through differentiation. Let us denote the standard basis of M"^ by 

ei = (l,0,---,0),e2-(0,l,0,...,0)---,e<i = (0,---,0,l). 

Theorem 1.9. Let d > 2. For a G Ng, n = \a\, define 

(1.9) p^ix) := J_lJ^||.||2|"l+<^-25" |||,||-<^+2| 

Then 

1. Pa E Tin (^^d Pa is the manic spherical harmonic of the form 

(1.10) pa{x)^x" + \\xfqa{x), qa E V^_2- 



2. Pa satisfies the recurrence relation 

(1-11) Pa+e,{x) = X,pa{x) - ^^ _ ^ H a^H ^5tPa(a^), 



3. {pa '■ \ot\ = n, ad = or 1} is a basis ofH'^. 

Proof. Taking the derivative of Pa{x) gives immediately the recurrence relation 
(1.11). Clearly po{x) =- 1. By induction, the recurrence relation shows that Pa is 
a homogeneous polynomial of degree n and it is of the form (1.10). We now show 
that Pa is a spherical harmonic. For g G Vf^ and p G R, a quick computation using 
Si=i Xidig{x) = ng{x) shows that 

(1.12) Adlxf.g) = p{2n + p + d- 2)\\x\\P-^g + \\x\\PAg. 

In particular, setting n = and g{x) — 1 gives A(||a;||^''+^) = 0. Furthermore, 
setting g = Pa and p = — 2n — d + 2 in (1.12) leads to 

^^"^^^ ^ 2"(d/2-l)J '^""°'^'''^°^^"^""'^'^ = °- 

Thus, Pq, e "H^. Since ||x|pq(x) is a linear combination of the monomials x^ with 
/3d > 2, by (1.10) and the hnear independence of {x" : \a\ = n,ad = orl}, 
it follows that the elements in the set {pa : la] = n, a^ = or 1} are linearly 
independent. The cardinality of the set is 

,. ^ri 1 ,. ^ri 1 fn + d~2\ fn + d—i 
d.mT^r + d.meH^ ,_2 JH d-2 

which is, by a simple identity of binomial coefficients and (1.5), precisely dirnHJ^. 
This completes the proof. D 

The right-hand side of (1.9) is called Maxwell's representation of harmonic poly- 
nomials ([8, 10]). The complete set of {pa ■ \oi\ — n} is necessarily linearly depen- 
dent by its cardinality. Moreover, by (1.9), 

r, -I- -l-T, — iziZl||-r-l|2|"l+rf-2/^aA /IItI|-'^+2"1 - 

Pa+2ei + ■ ■ ■+Pa+2ea — ——-JT-\\X\\ O n. \\\X\\ / — "Jj 

^ (,2 J" 

which gives dim7'„_2 linearly dependent relations among {pa ■ \a\ — n}. The set 
{pa '■ \ce\ — n} evidently contains many bases of "HJ^. The basis in item three of 
Theorem 1.9 is but one convenient choice. The proof of Theorem 1.9 relies on the 
fact that ||x|| ~'^+'^ is a harmonic function in R'* \ {0} for d > 2. In the case of d = 2, 
we need to replace this function by log||x||. Since the case d — 2 corresponds to 
the classical Fourier series, we leave the analogue of Theorem 1.9 for d = 2 to the 
interested reader. 

The basis {pa ■ \oi\ — n,ad ^ or 1} of H^ is not orthonormal. In fact, the 
elements of this basis are not mutually orthogonal. Orthonormal bases can be 
constructed by applying the Gram-Schmidt process. An explicit orthonormal basis 
for "H^ will be given in Section 5 in terms of spherical coordinates. 
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2. Projection operators and Zonal harmonics 

Let L^(S''^^) denote the space of square integrable functions on S''^^. Let 

proj„:L2(§rf-i)^7^^ 

denote the orthogonal projection from L^iW^^^) onto H'^. If P G Vf^ then P = 
Pn + \\xfQn, where F„ G Ui and 0„ e 'P^_2: by (1.4), so that proj„ P = P„. In 
particular, (1.11) shows that pa defined in (1.9) is the orthogonal projection of the 
function qa{x) — x"; that is, pa — proj„ qa- This leads to the following: 

Lemma 2.1. LetpeV^- Then 

(2.1) P-"J»^-E 4.,.(-n + 2-./2)," -"^^"^^"^- 

Proof. By linearity, it suffices to consider p being qa{x) — x" . By Theorem 1.9, 
proj„ qa{x) = Paix), and the proof amounts to showing that Pa{x) defined in (1.9) 
can be expanded as in (2.1). We use induction on n. The case n = is evident. 
Assume that (2.1) has been established for m = 0, 1, . . . , n. Applying (2.1) to qa(x), 
\a\ = n, it follows that 

a" {||a;r'^+2} =(-l)"2" (f - 1)^ ||a;||-2"-'^+2 

L"/2J 

j^^ 4^j!(-n + 2-rf/2)," " ^ 1 

Applying di to this identity, we obtain 

a^a" {||x|r''+2} = (-1)"2" (f - 1)^^ (-2n - rf + 2)||a;|r2"-'i+2 
L("+i)/2j , 

^ 4^i!(-n+l-d/2)," II ^ » t 1 7 ^^ n 

The terms in the square brackets are exactly A^jx^x"} and the constant in front 
simplifies to (_i)»+i2"+i (| - l)^^-^, so that the equation (2.1) holds for p{x) = 
Xix"-. This completes the induction. D 

Definition 2.2. The reproducing kernel Zn{-, ■) of Tif^ is uniquely determined by 

(2.2) — / Z„(x, y)p{y)da{y) = p(x), Vp e Ut xe E"-' 

and the requirement that Zn{x, •) be an element of Tif^ for each fixed x. 

That the kernel is well defined and unique follows from the Riesz representation 
theorem applied to the linear functional L{Y) :— Y{x), Y G H^, for a fixed x G 

Lemma 2.3. In terms of an orthonormal basis {Yj : 1 < J < dirnHJ^} ofH^, 

dim«^ 

(2.3) Z^{x,y)= Yl Yk(x)Yk{y), x,yeS''-\ 

and, despite (2.3), Z„ is independent of the particular choice of basis ofT-V^. 



Proof. Since Z„(x, •) G H'}^, it can be expressed as Zn{x,y) = '^k'^kYk{y) where 
the coefficients are determined by (2.2) as Ck = Yk(x). The uniqueness imphes that 
Zn is independent of the choice of basis. This can also be shown directly as follows. 
Let Y„ = (Yi, ...,Yn) with N — dimHf^ and regard it as a column vector. Then 
Zn{x, y) = [Y„(.x)]"'Y„(y). If {F/ : 1 < j < A^} is another orthonormal basis of Ut 
then Y^ — (5Y„. Since the orthonormality of {Yj} can be expressed as the fact that 
— Ld-i Y„(x)[Y„(x)]'''(i(7(a;) is an identity matrix, it follows readily that Q is an 
orthogonal matrix. Hence, Z^{x, y) = [Yn{x)Y'Q''QYn{y) - [Y'Jx)Y'Y'^{y). D 

The reproducing kernel is also the kernel for the projection operator. 
Lemma 2.4. The projection operator can be written as 

(2.4) proj„/(a;) = — / f{y)Z^{x,y)da{y). 

Proof. Since proj„ / G "HJ^, it can be expanded in terms of the orthonormal basis 
{Yj, 1 < j < Nn}, Nn ~ dimHfj, of H^, where the coefficients are determined by 
the orthonormality, 

Proj„/(a;) = Vcjy,(a;) with Cj = — f{y)Yj{y)da{y). 

~[ ^d Jgd-i 

If we pull out the integral in front of the sum, this is (2.4) by (2.3). D 

Lemma 2.5. The kernel Zn{-, •) satisfies the following properties: 

1. For every ^, r/ G E/^^^ , 

(2.5) — / Z„(e,y)Z„(77,y)da(y)=Z„(e,r7). 

^d Jgd-i 

2. Zn{x,y) depends only on (x, y). 

Proof. By Corollary 1.7, the uniqueness of Zn{x,y) shows that Zn{Qx,Qy) — 
Zn{x, y) for all Q e 0{d). Since for x,y £ S''"^ there exists a Q £ SO{d) such that 

Qx = (0, ..., 0, 1) and Qy = (0, ..., 0, a/1 — {x, y) , (x, y)), this shows that Zn{x, y) 
depends only on (x, y) . D 

From the second property of the lemma, Z„(x,y) — Fn{{x,y)), which is often 
called a zonal harmonic, since it is harmonic and depends only on {x,y). We now 
derive a closed formula for Fn , which turns out to be a multiple of the Gegenbauer 
polynomial, C^, of degree n defined, for A > and n E No, by 

(2.6) C„^(.):.(^,«,^,(-j;V;J^), 

where 2^1 is the hypergeometric function. The properties of the Gegenbauer poly- 
nomials are collected in Appendix B. 

Theorem 2.6. For n e No and x,y e S"^^^ , d>3, 

(2.7) Z„(x,y) = ^C^((x,y)), A = ^. 
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Proof. Let p e Hn- By Theorem 1.8, p{x) = (fc„(x, •),p)g. For fixed x, it follows 
from the same theorem that 

p{x) = (fc„(x, •),p)g = (proj„(fc„(x, •)),p)9 

= PT^oiJk,,{x,-)]{y)p{y)da{y). 

Since the kernel Zn{-^ ■) is uniquely determined by the reproducing property, this 
shows that Zn{x,y) = 2"(|)„ proj„[A:„(x, •)](?/). Since fc„(x, •) is a homogeneous 
polynomial of degree n and, taking the derivative on y, we have A^kn{x,y) — 
||x|p-'A:„_2j(2;, y), as is easily seen from 9ifc„(x, J/) — Xikn-i{x,y), Lemma 2.1 shows, 
for x,y e S''~\ that 

Z„(.x,y) = 2"(-) proj„[fc„(x,. )](?;)= J2 ^jf^ r-fc„-2,(x, y). 

Using the fact l/(n - 2j)! = {-n)2j/nl = 2'^i{-'j)j{^^)j/nl, we conclude then 

A n! ^ ]\{l-n-\)j 

= ^r^I^<^'^>^^^|^l-n-A=^' 

from which the stated result follows from (2.6). D 

Let {Yi : 1 < i < dimH^} be an orthonormal basis of "H^. Then (2.7) states 
that 

(2.8) Y. Y,{x)Y,{y)^^C^^{{x,y)), A = ^-. 

This identity is usually referred to as the addition formula of spherical harmonics, 
since for rf == 2 it is the addition formula of the cosine function (see Section 6). 

Corollary 2.7. For n e No and x,y e §'^"\ d>3, 

(2.9) \Znix,y)\<dimn'i and Z„(x,x) = dimH^. 

Proof. Set Fn{t) := ^C^,{t). By (2.7), Zn{x,x) ^ i^„(l) is a constant for all 
X € S'^"^. Setting x = y in (2.3) and integrating over S"^^^, we obtain 

1 /• 1 r ''""''" 

Fn{l) = — / Z„((a;,a;))dcr(x) = — / V if (2.)^^(2.) ^ dimH^. 

The inequality follows from applying the Cauchy-Schwarz inequality to (2.3). D 

Because of the relation (2.7), the Gcgcnbauer polynomials with A = ^^ are 
also called ultraspherical polynomials. A number of properties of the Gegenbauer 
polynomials can be obtained from the zonal spherical harmonics. For example, the 
corollary implies that C^(l) — ;jtj^ dimT^J^. Here is another example: 



Corollary 2.8. For A — ^^, the G eg enhauer polynomials C^ satisfy the orthogo- 
nality relation 

(2.10) ^^^ f C^{t)Ct{t){l ef-Ut = /.^5,„,„, 
where 

" n + X "^ ' 

Proof. Set again F„(i) = ^^^^C^{t). Since 2'„(x, •) and Zm{x, •) are orthogonal over 
S''""'^, and by (A. 5.1), their integrals can be written as an integral of one variable, 
we obtain 

^^^^ I F„(t)F„(t)(l - t^Y-^dt = — / F,,{{x,y))F„,{{x,y))da{y) 

^d J-i ^d Jgd-i 

- Fn{l)5rn,n = (dim H^)(5™,„, 

where the second line follows from (2.5) and Corollary 2.7. D 

The functions on S'^^^ that depend only on (x, y) are analogues of radial functions 
on R'' . For such functions, there is a Funk-Hcckc formula given below. 

Theorem 2.9. Let f be an integrable function such that J_^ |/(^)|(1 — t^Y'^^^^/'^dt 
is finite and d > 2. Then for every y„ G "H^, 

(2.11) / /((x,2/))y„(2/)da(y)=A„(/)y„(x), x e S'^-i, 
where Xnif) is a constant defined by 

Xnif) = Wd-l 



f\f{t)^iJ^{i-t-)'^dt. 

^-1 Cn' (1) 



■(1) 

Proof. If / is a polynomial of degree to, then we can expand / in terms of the 
Gegenbauer polynomials 

™ k + i^ d-2 

fc=0 2 

where A^ are determined by the orthogonality of the Gegenbauer polynomials, 

Afc = -^- f f{t)C^it)il - t^)'^dt, 
C\' (1)^-1 

and c^^ = /_]^(1 — t^)^^ dt = LOd/iOd-i- From (2.7) and the reproducing property 
of Zn(x, y), it follows that for n < m 

— f f{{x,y))Y,,{y)d<j{y) = A„y„(x), x e S^'K 

Since Xn/iOd ~ ^n{f) by definition, we have established the Funk-Hecke formula 
(2.11) for polynomials, and hence, by the Weierstrass theorem, for continuous func- 
tions, and the function satisfying the integrable condition in the statement can be 
approximated by a sequence of continuous functions. D 
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3. Zonal basis of spherical harmonics 

In view of the addition formula of spherical harmonics, one may ask whether 
there is a basis of spherical harmonics that consists entirely of zonal harmonics. 
This question is closely related to the problem of interpolation on the sphere. 

Throughout this subsection, we fix n, set A^ = dimH^ and fix {Yi, . . . , Fjv} as 
an orthonormal basis for "HJ^. Let {xi, . . . ,xjv} be a collection of points on S''^^. 
We let Ml :— li(xi) and, for fc = 2, 3, . . . , A^, define matrices 



Mk 



Yiixi) 



Ykixi) 



Yiixk)' 



Yk{xk)^ 



Mk{x) 





Yiix) - 


Mfc_i 






Yk^iix) 


|_ Yk{xi),...Yk{xu-i) 


Ykix), _ 



The product of M^ and its transpose Mj^ can be summed on applying the addition 
formula (2.8), as Mj^Mjy = [Zn{xi, Xj)]^^,^, which shows, in particular, 

(3.1) det[Z„(a;„Xj)]5^i = (detM^)^ > 0. 

This motivates the following definition. 

Definition 3.1. A collection of points {xi, 
system of degree n on the sphere S''^^ if 



. , X]\r} in § is called a fundamental 



det[C;J((x„x,))] 



N 



>o. 



A = 



d-2 
2 • 



Lemma 3.2. There exists a fundamental system of degree n on the sphere. 

Proof The existence of a fundamental system follows from the linear independence 
of {Yi, . . . , y^v}- Indeed, we can clearly choose xi E S'^^^ such that detMi = 
yi(a;i) 7^ 0. Assume that xi, . . . ,Xk, I < k < N — I, have been chosen such that 
det Mfe y^ 0. The determinant det Mk+i{x) is a polynomial of x and cannot be 
identically zero by the linear independence of {Fi, . . . , Yfe+i}, so that there is a 
Xk+i G S'^~^ such that det Mfe+i = det Mk{xk+i) y^ 0. In this way, we end up with 
a collection of points {xi , . . . , xn} on S'^~^ that satisfies det Mn ^ 0, which implies 
det [C;^((a;,,a;j))]f.^^ == CAr(det Ma,)^ > 0, where cat = A^/(n + A)^. D 



The proof shows, in fact, that there are infinitely many fundamental systems. 
Indeed, regarding xi,...,XAr as variables, we say that detM^r is a (rf — 1)A^- 
dimensional polynomial in these variables, and its zero set is an algebraic surface 
of M*^'^^^)^, which necessarily has measure zero. 



Theorem 3.3. 

then {Cl[{:,x,)) 

Proof Let P^{x) 



If{xi,... 
■■i = 1,2, 



,xn} is a fundamental system of points on the sphere, 

2 



,N},\ 



- , is a basis of Ti.^ 



^^C^{{-,x^)). The addition theorem (2 



8) gives 



P^{x) 



N 

E 
fc=i 



Yk{xi)Yk{x), 



1,2,. ..,Af, 



which shows that {Pi, . . . , Pn} is expressed in the basis {Yi, . . . , Iat} with transition 
matrix given by Mm = \Yk{xi)]^ j^^-^. Since {xi, . . . , xm} is fundamental, the matrix 
is invertible, by (3.1). We can then invert the system to express Yfc as a linear 
combination of Pi, ... , Pjv, which completes the proof. D 
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A word of caution is in order. The polynomial C^{{x,Xi)) is, for x G S''^^, a 
linear combination of the spherical harmonics according to the addition formula. 
It is not, however, a homogeneous polynomial of degree n in x G R''; rather, it is 
the restriction of the homogeneous polynomial ||x||"C;^((x/||x||, j/)) to the sphere. 
This is a situation in which the distinction between "HjJ and "H^^lgd-i is called for; 
see the discussion below Definition 1.1. 

Fundamental sets of points are closely related to the problem of interpolation. 
Indeed, it can be stated as follows: for a given set of data {{xj, yj) : 1 < j < N}, 
Xj e S''^^ and yj E R, there is a unique element Y E H^ such that Y{xj) = yj, 
j ~ 1, . . . ,N, if and only if the points {xi, . . . ,xn} form a fundamental system 
on the sphere. Much more interesting and challenging is the problem of choosing 
points in such a way that the resulting basis of zonal spherical harmonics has a 
relatively simple structure. 

A related result is a zonal basis for the space Vn of homogeneous polynomials of 
degree n and the space Ilf^ of all polynomials of degree at most n. 

Theorem 3.4. There exist points ^j^n G S''^^, 1 < .7 < ^"^ := dimP^, such that 
(i) {{x, £,j,n)" : 1 < j < r'l} is a basis for Vf^ of degree n. 

{a) For each polynomial f G 11^, there exist polynomials pj : [—1, 1] h- > R for 
^ 1^ j 1^ ffi such that 

fix) =Xl^j(^^'^J^"^)- 

Proof. Following the proof of the existence of fundamental system of points, it is 
easy to see that there exist points ^j,„ such that /(Cj,n) = if and only if / = for 
all / G Vn- For the proof of (i), we first deduce by the binomial formula that 

/,,„(x) := (x, 0,„)" - Yl S^^' "'^"' ^ - ^' - "■"• 

|q:|— n 

Let / e rl Then fix) = E|a|=„ a-^;". By (1.8), 

{fJi.n)a = "! XI '^"^j"" "" "■/(4»)' 

which implies, by the choice of ^j,„, that (/, fj,n)Q = 0, 1 < j < r^, if and only if 
/ = 0. Thus, {/j,„}-L = {0}. Since /j, „ G V^. This proves (i). 

For the proof of (ii), let m be an integer, < to < n — 1, and let fj.m^n '■— 
(x,^j,„)™. For / G V^, it follows as above that {f,fj,m.n)o = mlf{^j,m), 1 < j < 



r. 



Hence, if (/,/j,m,„)a = for 1 < j < <, then /(Cj,m) = (/5)te,m) = for 
1 < j < ^jj and all g G Vf^_.^, which implies by the choice of ^j^„ and the fact that 
Jg G V^, that fg = or f ^ 0. Consequently, V^ = span{/j,m,n ■ < j < r^} for 
< TO < n. Since 11^ — J2m=o'^m' this proves (ii). D 

4. Laplace-Beltrami operator 

The operator in the section heading is the spherical part of the Laplace operator, 
which we denote by Aq. The operator Aq plays an important role for analysis on 
the sphere. The usual approach to deriving this operator relies on an expression of 
the Laplace operator under a change of variables, which we describe first. 
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For X E K."^ let X 1-^ u = u{x) be a change of variables that is a bijection, so that 
we can also write x = x{u). Introduce the tensors 

Edxk dxk , i , s-^ dui duj , ^ ■ ■ ^ , 

fe=i ■' fe=i 

and let g := det{gi,j)f j^. Then (gij)^^ = (ff*'"')- A general result in Riemannian 
geometry, or a bit of tensor analysis, shows that the Laplace operator satisfies 

i=l « V y i^i j = i * 3 

The Laplace-Beltrami operator, i.e., the spherical part of the Laplace operator, can 
then be derived from (4.1) by the change of variables x i— >■ (r, ^i, . . . ,^d-i), where 
r > and ^ = (^i, . . . , ^d) G S'^^^. For this approach and a derivation of (4.1), see 
[10]. We shall adopt an approach that is elementary and self-contained. 

Lemma 4.1. In the spherical-polar coordinates x = r^, r > 0, ^ G S''^^, the 
Laplace operator satisfies 

d^ d-l d 1 



(4-2) A = — + — - + -Ao, 



where 

d—l i-j2 d—1 d—1 ^2 d—1 r. 

1=1 ^' j=i j=i ' ■' 1=1 ' 

Proof. Since ^ G S''^^, we have £,f + . . . + ^'^ = 1. We evaluate the Laplacian A 
under a change of variables {xi, . . . ,Xd) ^^ {f,^ii ■ ■ ■ jCd-i) under x = r^, which 
has inverse ^i — a;i/||a;||, . . . , S,d-i — Xd-i/\\x\\, r — \\x\\. The chain rule leads to 

oxi r ot,i r ^ o^j or 

(4-4) , , " 

_ ^d\^ f " Xd d 

dxd r^ ^ ■^ (9f o r dr ' 

j=i ■' 

If we apply the product rule for the partial derivative on x^, it follows that 

A=y7iA_livV- + e- ' 

,l( 1^'. 5 ^A" /. .2^^ 1 ^V 5 15 



where we have used Xd — ri^d, from which a straightforward, though tedious, com- 

c-2. 



putation, and simplification using ^^ + . . . +^^ — 1, establishes (4.2) and (4.3). D 



The Laplace-Beltrami operator also satisfies a recurrence relation that can be 
used to derive an explicit formula for Aq under a given coordinate system of S''"^. 
We write Ao,(i instead of Aq when we need to emphasize the dimension. 
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Lemma 4.2. Let Ao,d be the Laplacian-Beltrami operator for S''^^. For £, G S''^^, 
write ^ = (Vl -f^ri, t) with -1 < i < 1 and -q e S''"^. T/ien 

f-^) ^0. ^ ^jT^I (a - .')-!) + ^A.,_.. 

Proof. We work with the expression for Ao.d in (4.3) and make a change of variables 
{^i,---Ad-i) '-^ ijli, ■ ■ ■ ,Vd-2,t) defined by 



Ci = Vl - f^m^ ■■■, U-2 = Vl - i^?7d-2, Cd-i=t, 
where we have switched £,d-i and ^^ for convenience. The chain rule gives 

d Id , ^ d t ^ d d 

7^ = ~n=^^7^^ l<i<d-2, and = - r^Z^^j^ ^ J^^ 

d^i y/1 - 1^ orji o^d-i 1^^ "^ "Vj ot 

which can be used iteratively to compute Ao.d on writing (4.3) as 



d-2d-2 „2 <i-l Q 



A straightforward computation and another use of (4.3) then leads to 

a2 ,, , 9 1 

'W-^'-'^'m^—t^' 
which is precisely (4.5). D 



Ao,d = (1 - t')T^ - (rf - 1)^7^ + -, i^^o.d^i, 



The formula (4.3) gives an explicit expression for Aq in the local coordinates of 
S'^^^. An explicit formula for Aq in terms of spherical coordinates will be given in 
Section 5. 

Let V — (9i, . . . ,dd)- The proof of Lemma 4.1 also shows that 

(4.6) v-ivo+e|-, x = rC, ie^^-\ 

r or 

where Vq is the spherical gradient, which is the spherical part of V and involves 
only derivatives in ^. Its explicit expression can be read off from (4.4). We shall 
not need this expression and will be content with the following expression. 

Corollary 4.3. Let f e C^{^'^-^). Define F{y) ■- f{y/\\y\\), y G K'^. Then 

(4.7) Ao/(x)=AF(x) and Vof{x)^VF{x), x e S'^~\ 

The corollary follows immediately from (4.2) and (4.6), since a;/||x|| is indepen- 
dent of r. The expressions in (4.7) show that Aq and Vo are independent of the 
coordinates of S''^^. In fact, we could take (4.7) as the definition of Aq and Vq. 

The usual Laplacian A can be expressed in terms of the dot product of V, 
A = V • V, which can also be written-as is often in physics textbooks-as A = V^. 
The analogue of this identity also holds on the sphere. 

Lemma 4.4. The Laplace-Beltrami operator satisfies 

(4.8) Ao = Vo • Vo. 
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Proof. An application of (4.6) gives immediately 

1 1 / d\ d /I \ d^ 

A = V • V = — Vo • Vo + -Vo Kt^ + Ctt -Vo + jr^- 
j,z y, y (jr J or \r J or-^ 

We note that ^ • Vo/(0 — 0; since ^ e S''"^ is in the normal direction of ^ whereas 
Vo/(?)i t)y (4.6), is on the tangent plane at £_. Hence, we see that 

d fl \ 1 1 d 

^TT -^0 = --^ • Vo + -e • VoTT = 0. 
or \r J r^ r or 

Using (4.6), a quick computation gives Vo • ^ = d — 1, so that by the product rule, 

f d\ d d d 

Consequently, we conclude that 

92 d^l d 1 

^ = TT^ + IT + ^Vo • Vo. 

Comparing this with (4.2) completes the proof. D 

Our next result shows that the spherical harmonics are cigcnfunctions of the 
Laplace-Beltrami operator. 

Theorem 4.5. The spherical harmonics are eigenf unctions of Aq, 
(4.9) AoY{0 = -n{n + d-2)Y{0, ^Y eUt ^e^'^'K 

Proof. Let a; = r^, ^ e S''^^. Since Y G H^ is homogeneous, Y{x) — r"Y{£^) and 
by (4.2), 

= Ay (a;) == n{n - iy'-^Y{C) + {d - l)nr"-2y(^) + r'^-^AoYiO, 
which is (4.9) upon dividing by r"^^. D 

The identity (4.9) also implies that Aq is self-adjoint, which can also be proved 
directly and will be treated in the last section of this chapter, together with a 
number of other properties of the Laplace-Beltrami operator. 

5. Spherical harmonics in spherical coordinates 

The polar coordinates (xi,X2) = (r cos0,rsin0), r > 0, < 6 < 2tt, give 
coordinates for §^ when r = 1. The high-dimensional analogue is the spherical 
polar coordinates defined by 

xi — r sin 6d- 1 . . . sin 62 sin 9i , 
X2 ~ r sin 0rf_ 1 . . . sin 62 cos Oi , 
(5.1) 

Xd-i = r sin dd-i cos dd-2, 
Xd^rcosOd-i, 

where r > 0, < 0i < 27r, < ^i < tt for i = 2, . . . , d— 1. When r = 1 these are the 
coordinates for the unit sphere S''^^, and they are in fact defined recursively by 

X = {^smed-i,cos9d-i) e S'^-^, e e S'^-^. 

Let da — dad be Lebesgue measure on S''"^. Then it is easy to verify that 

(5.2) dadix) - ismed-if-^d9d-idad-iiO- 
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Since the Lebesgue measure of S^ is d9i, it follows by induction that 



d-2 

(5.3) da = d(Td = Y[ {sm9d-jf'^'^ dOd-i . . . d6'2d6'i 

in the spherical coordinates (5.1). Furthermore, (5.2) shows that 

(5.4) / f{x)d<Td{x)^ [ [ f{^ sine, cos 0)dad-i{£.){siney-^de. 
JS<i-i Jo Js<i-2 

The orthogonality (2.10) of the Gegenbauer polynomials can be written as 

(5.5) / C^(cos6')C^(cos6')(sin0)''-2^^^y — i_2_i;^A^^^^^ A = ^. 



r(|) 

Together with (5.4), this allows us to write down a basis of spherical harmonics in 
terms of the Gegenbauer polynomials in the spherical coordinates. 

Theorem 5.1. For d > 2 and a G Ng, define 

d-2 

(5.6) y„(x) := [/^„]-Vl"lg„(0i)[](sined_,)l"''''IC^^^(cos0rf_,), 

where ga{Oi) = cosad-iOi for ad — 0, sinad-i^*! for ad — l,\a^\ — aj + . . . + ad-i, 
Xj = |aJ+i| + {d-j- l)/2, and 

[rial — Oa J_J_ ._^ 

in which b^ = 2 if ad-i +ad > 0, while 6q, = 1 otherwise. Then {Yq : |a| — n,ad = 
0, 1} is an orthonormal basis of'H'^; that is, (Yq, Yg)gd-i — <^a.fj- 

Proof. To see that Y^ is a homogeneous polynomial, we use, by (5.1), the relation 
cos9k = Xk+i/ \ x\ + . . . + x"^^^ for 1 < /c < d — 1 to rewrite (5.6) as 

d-2 y 

„2 I , J2 \ai/2r^\i I Xd-j+1 



Y^{x) ^ [h^]-'g{x) Y{{xl + ... + x\_^^^r^l-'C% 



j = l ^^xl + . . . + X^_j^^ 

where g{x) = p"''-^ cosad-i0i for a^ = 0, p"''-'^ cos ad- lOi for ad = 1, with 
p = \/ x\ + x\. Since x\ — psin^i and X2 = pcos^i by (5.1), g{x) is either the real 
part or the imaginary part of {xi +ixi)"''~^, which shows that it is a homogeneous 
polynomial of degree ad-\ in x. Since C^{t) is even when n is even and odd when 
n is odd, we see that Y^ G T'^. Using (5.4), we see that 

(r„,r„.)g._i=^^^^-^ / ga{0i)9c.'{0i)dei 

^d Jo 



d-2 
X 



n / C^^,{cos9d-,)C^J{cos9d-,){smdd-,f^^d9d-j 



, = 1^0 



from which the orthogonality follows from the orthogonality of the Gegenbauer 
polynomials (5.4) and that of cosm9 and sinm^ on [0, 2tt), and the formula for ha 
follows from the normalization constant of the Gegenbauer polynomial. D 
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For d — 2 and the polar coordinates (xi,X2) — {rcos9,rsm9), it is easy to see 
that Vo — dg, where dg = d/d9. Hence by (4.8), the Laplace-Beltrami operator 
for d = 2 is Ag = dg. Using (4.5) iteratively, we see that the Laplace-Behrami 
operator Ag has an exphcit formula in the spherical coordinates (5.1), 



(5.7) Ao 



sm 



1 d 



Od-i dOd- 



d-2 



.d-2 r ^ 



sm tld-i- 



yd-i 



d-2 

E- 



■f sin 9d-i ■ ■ ■ sin 9j+i sin-' 9j d9j 



'09, 



6. Spherical harmonics in two and three variables 

Since spherical harmonics in two and three variables are used most often in 
applications, we state their properties in this section. 

6.1. Spherical harmonics in two variables. For d — 2, dimH^ — 2. An or- 
thogonal basis of "H^ is given by the real and imaginary parts of (xi + ix2)^ , since 
both are homogeneous of degree n and are harmonic as the real and imaginary 
parts of an analytic function. In polar coordinates (a;i,a;2) — {r cos9,rs,va9) of M^, 
this basis is given by 

(6.1) YnHx) = r"cosn6l, YJi^^x) = r" sin 716*. 

Hence, restricting to the circle S^, the spherical harmonics are precisely the cosine 
and sine functions. In particular, spherical harmonic expansions on S^ are the 
classical Fourier expansions in cosine and sine functions. 

As homogeneous polynomials, the basis (6.1) is given explicitly in terms of the 
Chebyshev polynomials r„ and t/„ defined by 

sin(n + \)9 
TriU) — cosn9 and Unit) = , where t — cos 9, 

which are related to the Gegenbauer polynomials: t/„(t) = Cni^) ^^d 

lim ic;J(x) = -T„(x). 
A->-o+ A n 

The basis in (6.1) can be rewritten then as 

(6.2) yW(x)=r"r„(^), Y('Hx)^r--'x,U^-i{y), 



which shows explicitly that these are homogeneous polynomials since r = \/ x\ + x\ 
and both T„(i) and C/„(t) are even if n is even, odd if n is odd. 

When d — 2, the zonal polynomial is given by Tn{{x, y)) — cosn{9 — 4>) and the 
addition formula (2.8) becomes, by (6.1), the addition formula 

cosn9 cosncf) + smn9smn(j) — cosn{9 — (f)). 

The expression (4.2) of the Laplace operator in polar coordinates becomes 

A-— -— — — ■ 
dr r dr r^ d9'^ 

in particular, the Laplace-Beltrami operator on §^ is simply Aq = d^ /d9'^. 
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6.2. Spherical harmonics in three variables. The space H^ of spherical har- 
monics of degree n has dimension 2n + 1 . For d = 3 the spherical polar coordinates 
(5.1) are written as 

xi — r sin^sini/), 

(6.3) { X2 ^r sine cos (j), , < 6* < tt, < < 27r, r > 0. 
x^ — r cos 6, 

The surface area of S^ is Air and the integral over S''^^ is parameterized by 

(6.4) / f{x)da= / / /(sin6'sin(/), sin0cos(/), cos0)(i(/)sin0(i0. 

JS2 Jo Jo 

The orthogonal basis (5.6) in spherical coordinates becomes 

Yki{d,(f>) = (sin 0)'=C^+| (cos 6*) cos fc0, < k < n, 

(6.5) ' J 

^fe"2(^,'/') == (sin 0)'=C^+| (cos 6*) sin fc0, 1 < k < n. 

Their i^(§^) norms can be deduced from (5.6). This basis is often written in terms 
of the associated Legendre polynomials P^ (t) defined by 

P,':ix) := (-1)"(1 - x'f/'-^P4x) = (2fc - 1)!!(-1)"(1 - x'f/'c'^t'J'ix), 

where P„ (t) ~ C° (t) denotes the Legendre polynomial of degree n (see Appendix B 
for properties of P„ and P^), and in terms of {e*'^'^, e^^'"^} instead of {cos fc0, sinfc^}. 
In this way an orthonormal basis of 'H'^ is given by 

(6.6) Y,M9,^) = ( ^^"|„+^fc)7^^T ^'^"(""^^)"''"' -^^fc^"- 

The addition formula (2.8) then reads, assuming that x and y have spherical coor- 
dinates {0,(f>) and [9' ,(p'), respectively, 

n 

(6.7) Y. YkAd,(l>)Yk,n{e'A') = {'2n+l)Pn{{x,y)). 

k—~n 

In terms of the coordinates (6.3), the Laplace-Beltrami operator is given by 

Id/ d\ 1 92 






sin^ e d(j>' 
as seen from (5.7). 

7. Representation of the rotation group 

In this section we show that the representation of the group SO{d) in spaces of 
harmonic polynomials is irreducible. 

A representation of SO{d) is a homomorphism from SO{d) to the group of 
nonsingular continuous linear transformations of L'^(S'^^^). We associate with each 
element Q E SO{d) an operator T(Q) in the space of L'^{E>'^^^), defined by 

(7.1) T{Q)f{x) - f{Q-^x), x e S-^-i. 

Evidently, for each Q e SO{d), T{Q) is a nonsingular linear transformation of 
L'^{§'^^^) and T is a homomorphism, 

T{QiQ2) - TiQi)T{Q2), VQi, O2 e SO{d). 
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Thus, T is a representation of SO{d). Since da is invariant under rotations, 
l|r(<3)/||2 = II/II2 in the L2(§'^-1) norm, so that T{Q) is unitary. 

A Hnear space U is cahed invariant under T if T{Q) maps U to itself for all 
Q e SO{d). The null space and L^(§''~^) itself are trivial invariant subspaces. A 
representation T is called irreducible if it has only trivial invariant subspaces. The 
space "HJ^ is an invariant subspace of T in (7.1). 

Let Tnd denote the representation of SO{d) corresponding to T in the invariant 
subspace "HJ^. We want to show that Tn^d is irreducible. 

Lemma 7.1. A spherical harmonic Y G Ti^ is invariant under all rotations in 
SO{d) that leaves Xd fixed if and only if 

(7.2) r(x)=c||x|rC„^(^), A=^. 

where c is a constant. 

Proof. If Y is invariant under rotations that fix Xd and is a homogeneous polynomial 
of degree n, then it can be written as 

0<j<n/2 0<j<n/2 

where the second equal sign follows from expanding (||a;|p — x'^Y and changing the 
order of summation. Since Y is harmonic, computing AY{x) = shows that Cj 
satisfies the recurrence relation 

4(j + l)(n - J - l)cj+i + (n - 2j){n - 2] - l)c, = 0. 

Solving the recurrence equation for c,- we conclude that 



r(x) = co E 



y 2)3\ 2 I] »-2j|i ||2.7 



x^ -lla;!! 



Consequently, (7.2) follows from the formula (??) for Gegenbauer polynomials. 
Since the function Y in (7.2) is clearly invariant under all rotations that fix Xd and 
we have just shown that it is harmonic, the proof is complete. D 

Theorem 7.2. The representation Tn,d of SO{d) on "H^ is irreducible. 

Proof. Assume that U is an invariant subspace of H^ and U is not a null space. 
Let {Yj : 1 < j < M}, M < dimH^, be an orthonormal basis of U. Follow- 
ing the proof of Lemma 2.5, there is a polynomial F{t) of one variable such that 
J2j=i '^j{^)^j{y) — ^{{^^ y))- Ii^ particular, setting y — Cd — (0, ..., 0, 1) shows that 
F{{x, Cd)) is in "H^ and it is evidently invariant under rotations in SO{d) that fix Xd. 
Hence, by Lemma 7.1 F{{x,ed)) = c\\x\\^C^{j4t). In particular, this shows that 

||x||"C^(^^) G U. On the other hand, let U-^ denote the orthogonal complement of 

U in nt If / e U^ and g G U, then (T(0)/, g)g._i = {f,T{Q-^)g)s,-i = 0, which 
shows that U^ is also an invariant subspace of H"^. Applying the same argument 
as for W shows then \\x\\"^C'^{^A €U^, which contradicts to Z^ n Z^-^ = {0}. Thus, 
U must be trivial. D 
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8. Angular derivatives and the Laplace-Beltrami operator 

Consider the case d — 2 and the polar coordinates (a;i,X2) = {r cos 9, r sin 9). 
Let dr = d/dr and de = d/d9, while we retain di and 82 for the partial derivatives 
with respect to xi and X2. It then follows that 

sin 9 ens f) 

di = cos9dr dg, 82 = sm9dr 



r r 

From these relations it follows easily that the angular derivative dg can also be 
written as dg = a;i92 — a:;2C^i a-nd the operator Aq is Aq — dg. We introduce angular 
derivatives in higher dimensions as follows. 

Definition 8.1. For x eM."^ and 1 < i y^ j < d, define 

d 

(8.1) Aj := x,dj - Xjd^ = ^^, 

where 9ij is the angle of polar coordinates in the {xi,Xj) -plane, defined by {xi,Xj) — 
rij{cos6i,j,sui9ij), ri,j > and < 9ij < 2tt. 

By its definition with partial derivatives on R'', Dij acts on M"^, yet the second 
equality in (8.1) shows that it acts on the sphere S''"^. Thus, for / defined on R'*, 

(8.2) (A,,/)(0 = A,,[/(0], eeS-^-i, 

where the right-hand side means that Dij is acting on /(^). 

Since Dji = —Dij, the number of distinct operators Dij is (2). The operator 
Aq can be decomposed in terms of them. 

Theorem 8.2. On S''^^, Aq satisfies the decomposition 

(8.3) ^0= J2 ^Ir 

l<i<j<d 

Proof. Let F{x) — f{ TTfir ) • A straightforward computation shows that 

d^F d - 1 ^ 
dxidxj \\x\\ ^-^ dx 



l<i<j<d " " i=l j = l 



Consequently, restricting to § and comparing with (4.3), (8.3) follows. D 

Let Qi,j,g denote a rotation by the angle 9 in the (a;i,Xj) -plane, oriented so 
that (xi,Xj) = (scos9, ssm9). Then T{Qi_j^g), defined in (7.1), maps / into 
T{Qij,g)f{x) — f{Qij^^gx). Written exphcitly, for example, for (i,j) = (1,2), 
we have 

(8.4) T{Qi^2.9f){x) — f{xi cos 9 + X2 sin0, —xi sin 9 + X2 cos 9, X3, . . . , Xd). 

Then Di^ is the infinitesimal operator of T{Q), 



/g 5\ dT{Qij^g) 



d9 



d A- n 

?=0 dXn OXi 



"^0 

where the first equality follows from (8.4). The infinitesimal operator plays an 
important role in representation theory; see, for example, [12]. 

The operators Dij will play an important role for approximation theory on the 
sphere. We state several more properties of these operators. 
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Lemma 8.3. For 1 < i < j < d, the operators Dij commute with Aq. In particu- 
lar, Dij maps Tin to itself. 

Proof. By symmetry, we only to show only that I?i.2- Let [^, -B] = AB — BA 
denote the commutator of A and B. A quick computation shows that 

[dj,Dk,i] = 5j,kdi - Sj.idk, [xi, Dk.i] = Si^kXi - Si^iXk, 
from which it is easy to see that, for 1 < i < j < d, 1 < k < I < d, we have 

(8.6) [Aj, -Dfe,;] = -Si^kDjj + SijDj^k + ^j.kDi^i - SjjDi^k- 

Using (8.6), a simple computation shows that [Di,2,Di i] = —(D1.1D2.1 + D2,iDi.i) 
and [Z?i,2,^2,;] = Di.iD2,i + -D2,;£'i,/, so that [^,2,1^?,/ + ^Ll = for i > 2. 
Moreover, by (8.6), [Di2,Df i] — whenever 3 < k < I < d. Summing over {k,l) 
ioT 1 < k < I < d then proves [£'1,2, Aq] = 0. D 

Proposition 8.4. For f,g e C''^{S'^-'^) and 1 < i ^ j < d, 

(8.7) / f{x)D,,jg{x)da{x) - - / D,,jf{x)g{x)d(j{x). 

Proof. By the rotation invariancc of the Lebesgue measure da, we obtain, for every 
0€ [-TT,Tr], 

f{x)g{Qi.j-ex)da{x)= / f{Q,,j^0x)g{x)da{x). 

Differentiating both sides of this identity with respect to 6 and evaluating the 
resulting equation at 6* — leads to, by (8.5), the desired equation (8.7). D 

The equation (8.7) allows us to define distributional derivatives _D[ on S''^^ for 
r € N via the identity, wtih g G C°°(§''~i), 

Dljix)gix)daix) = {-ly f f{x)Dl^g{x)da{x). 

Summing (8.8) with r = 2 over i < j and applying (8.3), it follows immediately 
that Aq is self-adjoint, which can also be deduced from Theorem 4.5. 

Corollary 8.5. For f,g e C^{^'^-'^), 

f{x)Aog{x)d(T = / Aof{x)g{x)da. 

The spherical gradient Vq is a vector of first order differential operator on the 
sphere, which can be written in terms of Dij as follows. 

Lemma 8.6. For f E C^(§'^^^) and I < j < d, the j-th component 0/ Vq/ satisfies 

(8.9) (Vo),/(0= E ?^^^j-/(0, ee§'-^ 

l<i<d,i^j 

Furthermore, for f,g G C^{E>'^^^), the following identity holds: 

(8.10) Vo/(e) • Vo5(C) = E Aj/(0A,,5(e), ieS"-'. 

l<i<j<d 



5.11) (Vo),/(0 = j^FiO = d,f - c, Y. C.9./ 
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Proof. Let F(j/) = /(2//||y||)- From the definition and ||^|| = 1, we obtain 

d 

dxj 

«— 1 2—1 -i— 1 

which gives (8.9) since Di^if — 0. (8.11) means that 
(8.12) Vo/(0 = V/(0-e(e-V/). 

Since ^ • Vo/(0 = 0^ using (8.11) and the definition of -Di.j, it fohows that 

d d d 

Vo/ • Voff - Vo/ • V5 - E ^J/^J3 " ^ ^ ^^iAf^l9 = E A,j/(0A,,5(0, 

j = l j = l i=l i<j 

where the second equahty uses ||^|| = 1. D 

As an apphcation, we state an integration by parts formula on the sphere. 
Proposition 8.7. For f,g e C^{S'^-^), 



(8.13) / .f{x)Vog{x)da = - / (Vo/(a;) - {d - l)xf{x)) g{x)da. 
Furthermore, for J G C^{^'^-^) and g G C\S'^-^), 

(8.14) / Vo/ • Vogda = - f Aof{x)g{x)da. 
Proof. Using the expression (8.11) and applying (8.7), we obtain 

f{x){\/n)jg{x)da = E / Xif{x)Di,jg{x)da 

1 J^d-l 

l<i<d,iyij''^ 

= - E / D^^j{x,f{x))g{x)da. 

By the chain rule, recalling (8.2) if necessary, we have Dij{xif{x)) — XiDijf — 
Xjf{x). Hence we obtain 

/ f{x){S/o)0{x)da = - f ( E ^^D^^,f{x)-{d-l)f{x))g{x)da, 

which gives the j-th component of (8.13). Equation (8.14) follows immediately 
from (8.10), (8.7) and (8.3). D 

9. Notes 

Spherical harmonics appear in many disciplines and in many different branches of 
mathematics. Many books contain parts of the theory of spherical harmonics. Our 
treatment covers what is needed for harmonic analysis and approximation theory 
in this book. Below we comment on some books that we have consulted. 

A classical treatise on spherical harmonics is [8], a good source for classical 
results. A short but nice expository is [9], which was later expanded into [10]. The 
reference book [5] contains a chapter on spherical harmonics. A rich resource for 
spherical harmonics in Fourier analysis is [11]. Applications to and connections 
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with group representation arc studied extensively in [12]; see also [7]. For their 
role in the context of orthogonal polynomials of several variables, see [4] as well 
as [5]. The book [1] contains a chapter on spherical harmonics in light of special 
functions. The theory of harmonic functions is treated in [3], including material 
on spherical harmonics. The book [6] deals with spherical harmonics in geometric 
applications. Finally, the recent book [2] provides an introduction to spherical 
harmonics and approximation on the sphere from the perspective of applications in 
numerical analysis. 

Aside from their role in representation theory, the operators Dij do not seem to 
have received much attention in analysis. Most of the materials in Section 8 have 
not previously appeared in books. These operators play an important role in our 
development of approximation theory on the sphere. 
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